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INTRODUCTION 


This  report  is  a  tutorial  for  readers  interested  in  the  effects  of 
vibration  on  the  output  frequency  of  a  quartz  crystal  controlled  oscillator. 
It  is  a  compendium  of  information  gathered  from  various  published  sources 
plus  several  experimental  examples.  Individuals  seekinf  more  details  are 
encouraged  to  refer  to  the  works  cited  in  the  references. 

The  acceleration  sensitivity  of  crystal  oscillators  is  assuming 
increasing  significance  in  modern  communications,  navigation  and  identi¬ 
fication  systems.  The  requirements  that  these  systems  be  secure  and 
resistant  to  jamming  imposes  stringent  requirements  on  the  frequency  control 
aspects  of  the  system  designs.  For  example,  in  systems  such  as  the 
NAVSTAR  Global  Positioning  System  (GPS),  the  vibration  sensitivity  of  the 
crystal  oscillator  is  a  major  limitation  on  the  jam  resistance  achievable 
in  a  mobile  (i.e.,  vibratory)  environment.1’2 

The  resonant  frequency  of  a  quartz  crystal  resonator  is  influenced 
by  acceleration. 3  A  classic  example  of  this  is  the  static  "2g  tip-over" 
test,  whereby  a  resonator  exhibits  a  variation  in  resonant  frequency 
depending  on  its  orientation  with  respect  to  the  earth's  gravitational 
field.  The  rotation  of  a  resonator  by  180°  about  an  axis  parallel  to  the 
surface  of  the  earth  causes  a  reversal  of  sign  and  therefore  a  change  in 
the  gravitational  acceleration  vector  of  amplitude  2g  (g  being  the 
acceleration  due  to  gravity  at  sea  level).  The  source  of  the  coupling 
between  the  frequency  and  the  gravitational  field  is  most  probably  the 
stress  sensitivity  of  quartz.4  A  stress  applied  to  a  quartz  plate  will 
change  the  resonant  frequency  of  that  plate.  The  effect  is  linear  for 
small  stresses.5  A  variation  in  the  acceleration  direction  will  vary  the 
stress  distribution  and  therefore  the  frequency.  The  "2g  tip-over"  test 
demonstrates  that  a  change  in  sign  of  the  acceleration  vector  with  respect 
to  the  resonator  will  change  the  sign  of  the  frequency  shift. 

The  acceleration  due  to  gravity  (even  of  a  static  resonator)  is 
indistinguishable  from  an  inertial  acceleration. 6  Therefore,  a  resonator 
undergoing  inertial  acceleration  will  exhibit  a  frequency  shift  with  a 
magnitude  dependent  on  the  acceleration  magnitude  and  sign  dependent  on 
the  direction  of  the  acceleration. 

SIMPLE  HARMONIC  MOTION 


Consider  a  resonator  (presumably  in  some  sort  of  oscillator  circuit 
so  that  one  can  monitor  the  resonant  frequency)  undergoing  simple  harmonic 
motion  (i.e.,  a  simple  one  frequency  vibration).  The  position  of  the 
resonator  as  a  function  of  time  may  be  written  as 


X(t)  =  (AX  / 2)  cos ( 2tt f  t) 
max  v 


(1) 


X ( t )  =  position  as  a  function  of  time 
1 


where 


■ 


f  =  vibration  frequency 
^max  =  peak-to-peak  displacement. 

Let  A(t)  be  the  acceleration  as  a  function  of  time.  Since 

A(t)  =  d2X(t)/dt2  (2) 


we  get 


A(t)  =  ~ ( 2tt f  t)2  (AX  12)  cos  (2irf  t)  (3) 

v  max  v 


or 


A(t)  =  -(A  )  cos  ( 27rf  t) 

max  v 


(3a) 


where 


A 


max 


(2trf  t)2  (AX  12) 
v  max 


(3b) 


If  a  resonator  experiences  such  a  motion,  the  frequency,  F0,  of  the  output 
signal  will  be  modulated  in  phase  with  the  acceleration,  due  to  the  time 
varying  stress,  by  an  amount  proportional  to  the  acceleration  magnitude. 
The  frequency  excursion,  AF(t),  is  then 

AF(t)  «  A(t)FQ  (4) 


where  the  constant  of  proportionality  is  known  as  the  vibration  sensitivity, 
which  will  be  denoted  by  y.  The  maximum  frequency  excursion  is  then  related 
to  the  acceleration  by  the  following  relation: 


AF  =  A  yF  . 
max  max  o 


(5) 


The  output  frequency  is, therefore. 


F(t)  =  F  +  AF  cos(27rf  t) 
o  max  v 


(6) 


2 


where 


F(t)  =  output  frequency 


Fq  =  "at  rest"  frequency 


Figure  1  is  a  representation  of  the  output  frequency,  F(t),  at 
different  times.  It  can  be  seen  that  the  instantaneous  value  of  F(t)  is 
modulated  by  an  amount  AF  about  F  at  a  rate  of  f  .  The  system  designer 
must  take  this  modulation  into  account.  A  narrow  bandwidth  servo-loop 
can  lose  lock  if  the  frequency  deviates  too  far  from  center. 

The  output  voltage  of  the  monitoring  oscillator  is 


V(t)  =  V  cos  $(t) 
o 


(7) 


where  <J>(t),  the  "phase",  is  related  to  the  frequency  of  the  resonator  by 


t 

4>(t)  =  2tt  /  F(t)  dt  (8) 

o 


therefore 


t  t 

$(t)  =  2tt  f  F  dt  +  2tt  f  AF  cos ( 2 tt f  t)  dt  (9) 

4  o  4  max  v 


and 


$(t)  =  2ttF  t  +  (AF  /f  )  sin(2rrf  t)  (10) 

o  max  v  v 


whence  it  follows 


V(t)  =  V  cos  (2uF  t  +  (AF  /f  )  sin(27rf  t))  (11) 

o  o  max  v  v 


Figure  2  is  an  example  of  the  output  voltage  of  a  frequency  mod¬ 
ulated  resonator.  The  low  frequency  curve  is  the  modulating  frequency 
(vibration) . 

It  is  very  difficult  to  quantify  data  in  the  form  displayed  in 
Figure  2.  This  representation  is  known  as  the  time  domain  representation 
(the  horizontal  axis  is  time).  A  much  clearer  picture  can  be  gained  by 


TIME 


presenting  the  same  information  in  the  frequency  domain  where  the  output 
voltage  of  the  oscillator  is  character ized  as  a  function  of  frequency. 

To  measure  the  output  of  an  oscillator  in  the  frequency  domain,  one  must 
measure  the  signal  passed  through  a  series  of  narrowband  transmission 
filters.  One  can  perform  this  task  most  conveniently  with  a  spectrum 
analyzer.  In  mathematical  terms  the  frequency  domain  output  of  the 
spectrum  analyzer  corresponds  to  the  decomposition  of  a  complex  waveform 
into  its  frequency  (Fourier)  components. 

FREQUENCY  DOMAIN 


The  decomposition  of  V(t)  into  its  frequency  components  is  most 
simply  performed  with  the  aid  of  several  trigonometric  identities.7 


cos(X  +  Y)  -  cos(X)  cos(Y) 

cos(X)  cos(Y)  =  'i  lcos(X  +  Y) 

-sin (X)  sin ( Y)  =  tcos(X  +  Y) 

and  two  relations  derived  in  Appendix  I 

cos(Ssin  X)  =  J  (6)  +2  j 
o  i 

n=l 

OO 

sin (Bsin  X)  =  2  £  J ^  + 


n=o 

where  Jp($)  =  Bessel  function  of  order  p 
From  equation  11  above 


-  sin(X)  sin(Y) 

+  cos (X  -  Y)} 

-  cos (X  -  Y)} 

( 3 )  cos(2nx) 

(B)  sin((2n  +  l)x) 
(see  Appendix  1)  and 


(12) 

(13) 

(14) 


(15) 

(16) 


argument  B. 


V(t)  =  V  {cos  (2ttF  t  +  (AF  /f  )  sin(2Trf  t))} 
o  o  max  v  v 

using  equation  12  with  X  =  2ttF  t  and  Y  =  (AF  /f  )  sin(2TTf  t) 

.  o  max  v  v 

we  get 

V(t)  =  V  (cos(2ttF  t)  cos  ( (AF  /f  )  sin(2iTf  t)) 
o  o  max  v  v 

-sin(2irF  t)  sin((AF  /f  )  sin(2irf  t))} 
o  max  v  v 

using  equations  15  and  16  with  B  =  (AF  /f  )  and  X  =  2 it f  t  we  get 

max  v  v 


(17) 

(18) 


(19) 
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V(t)  =  V  {cos  (2tiF  t)(J  (AF  /  f  ) 
o  o  o  max  v 

oo 

+  2  £  J„  (AF  /f  )  cos(4niTf  t)) 

,  2n  max  v  v 

n=l 

OO 

-sin(2uFot)(2  l  J2n+i(AFmax/fv-)  sin ( (2n+l )  2irfvt))}  (20) 
n=o 

which  reduces  to 


V(t)  =  V  {J  (AF  / f  )  cos (2nF  t) 
o  o  max  v  o 


+  2  [  J,  (AF  /f  )  cos(2nl  t)  cos(4irnf  t) 

L,  2n  max  v  o  v 

n=l 


"  2  E  J2n+l(AFmax/fv)  sin(2TrF  t)  sin((2n+l)  2irf  t)} 
n=o 


(21  ) 


using  equation  13  with  X  =  2nF  t  and  Y  =  4nirf  t  and  equation  14  with 
X  =  2ttF  t  and  Y  =  (2n+l)  2rr f  ?  we  get 


Y(t)  =  V  {J  (AF  /f  )  cos  (2ttF  t) 
o  o  max  v  o 


+  ^  J2n^AF™o„/f„^cos(27T(F^+2nf„)t)  +  cos  (2tt (F^-2nf<;)  t) 


n=l 


max  v 


+  l  J2n+l(AFiaax/fv)(cos(2Tr(Fot(2n+l)f  )t) 
n=o 


-  cos  (2tt  (F  -  (2n+l)  )  f  t } 


(22) 


which  expands  to 


V(t)  =  V  {J  (AF  /  f  )  cos  (  2tt  f  t)  +  J,  (AF  /  f  )  cos(2tt(F  +f  )t) 
o  o  max  v  v  1  max  v  o  v 

-  J.(AF  /  f  )  cos  (2m(F  -f  )t)  +  J0(AF/f  )  cos(2-(F  +2  f  )t) 

1  max  v  ov  2  v  o  v 

+  J0(AF  /f  )  cos (2i(F  -2 f  )t)  +  .  .  .  .  r  (23 

2  max  v  o  v 

So  the  spectrum  of  a  sinusoidally  vibrated  resonator  contains  an 
infinite  number  of  Fourier  components  on  each  side  of  the  carrier, 
separated  from  the  carrier  by  f^,  2fv>  3fv  etc.  Each  of  these  terms 
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is  called  a  vibration  induced  sideband. 

Example  1: 

The  acceleration  sensitivity  of  AT-cut  quartz  is  typically  about 
2  x  10  9/g.  For  a  resonator  with  a  frequency  of  5  MHz,  undergoing  a 
harmonic  vibration  with  maximum  acceleration  of  5g,  the  AF  is 


AF  =  yA  F  =  (2  x  10~9/g) (5g) (5  x  106  Hz)  (24) 

max  max  o  a  b 

AF  =  0.05  Hz 
max 


If  the  vibration  frequency  is  above  5  Hz, 
is  below  0.01.  As  can  be  seen  in  Appendix  II, 
<0.1. 


the  value  of  AF  /f 

for  values  of  A?ax  /f 
max  Lv 


J  (AF  /f  <0.1)  -  1 
o  max  v 


J  (AF  /f  <0.1)  =  (AF  /f  )/2 
o  max  v  max  v 


J  .  . (AF  /f  <0.1)  -  0 
P>1  max  v 


(25) 


Therefore,  only  the  first  upper  and  lower  sidebands  are  important  for 
small  modulations.  Mathematically,  this  reduces  equation  23  to 


V(t)  -  V  {J  (AF  /f  )  cos (2ttF  t)  +  J,  (AF  /f  )  cos(2tt(F  +f  )t) 
o  o  max  v  o  1  max  v  o  v 


-  J. (  F  / f  )  cos (2  (F  -f  )t)} 
1  max  v  o  v 


(26) 


let  n(fv)  be  defined  as  the  ratio  of  the  power  contained  in  the 
n  vibration  induced  sideband  to  the  carrier  power.  That  is 

V  )  =  (J.  (AF  /f  )/J  (AF  /f  )}2  (27) 

^  v  v  1  max  v  o  max  v 

to  express  in  decibels 


or 


izf  v“(fv)(dB)  -  10  log  <Sfv"(fv> 

^v"(fv)(dB)  -  20  log{JlWFmax/fv)/Jo«Fmji/fv)) 


(28) 

(29) 
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For  AF  /f  <0.1 
max  v 


c^fvn(fv)(dB)  =  20  log  (AFmax/2fv) 


(30) 


For  our  current  example: 

o2fv1  (1  Hz)  =  -  32dB 

o2f  1  (5  Hz)  =  -  46dB 

V 

(25  Hz)  =  -  60dB 

(50  Hz)  =  -  66dB 

c^1  (500  Hz)  =  -  86dB 

c^f  1  (5000  Hz)  =  -106dB 

The  negative  sign  means  that  the  sideband  is  smaller  than  the 
carrier.  In  the  25  Hz  case,  for  example,  the  sideband  is  smaller  by  a 
factor  of  1  x  106  in  power,  or  a  1  volt  carrier  signal  will  have  a  1 
millivolt  sideband. 

It  is  seen  that  the  lower  the  vibration  frequency  the  smaller  the 
difference  between  the  carrier  and  the  sideband  amplitudes,  i.e.,  the 
signal  to  noise  ratio  is  degraded.  Unfortunately,  the  low  frequency 
vibrations  are  the  ones  most  difficult  to  damp. 

Figure  3  is  an  actual  trace  of  the  power  spectrum  of  a  resonator  at 
rest.  The  zero  frequency  point  on  the  abcissa  is  the  resonant  frequency, 
F  .  The  abcissa  is  the  Fourier  frequency  minus  the  oscillator  "at  rest" 
frequency . 

Figure  4  is  the  power  spectrum  of  the  same  resonator  being  vibrated 
at  a  level  of  lOg  at  a  vibration  frequency  of  100  Hz.  The  lobes  on 
either  side  of  the  carrier  are  the  vibration  induced  sidebands.  The 
ordinate  is  in  units  of  dB  below  the  carrier.  It  can  be  seen  that  the 
ratio  of  the  power  in  the  sideband  to  that  of  the  carrier  isj-60.1  dB. 
This  translates  into  an  acceleration  sensitivity  of  4.0  x  10  Vg. 

FREQUENCY  MULTIPLICATION 

In  systems  use,  it  is  common  for  an  oscillator  output  frequency 
to  be  multiplied  to  some  higher  frequency.  Consider  a  vibration 
modulated  frequency,  as  in  the  previous  section. 
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RES  BW  10  H* 


STATIC'  Rl 


FIGURE  4  -  POWER  SPECTRUM  OF  A  RESONATOR  UNDERGOING 
108Hz.  10g  VIBRATION 


(32) 


F(t)  =  F  +  AF  cos(2nf  t) 
o  max  v 


If  this  frequency  is  multiplied  by  a  factor,  N,  the  new  frequency, 

F  (t),  is  then 
m 


F  (t)  =  NF(t)  =  N(F  +  AF  cos(2irf  t)) 
m  o  max  v 

As  earlier,  the  phase  ( t )  is  given  by 

t  t 


<t>  (t)  =  2tt  /  F  ( t)dt  »  2ttN  /  F(t)dt 

m  o  m  o 

We  have  also  seen  that  this  may  be  written  as 

<J>  (t)  =  2ttNF  t  +  (NAF  /f  )  sin(27rf  t) 
m  o  max  v  v 

so  V ( t )  =  V  cos(2^NF  t  +  (NAF  /f  )  sin(2nf  t)), 

o  o  max  v  v 

From  this  it  follows  that 


(33) 


(34) 


(35) 

(36) 


V(t)  *  V  {J  (NAF  /f  )  cos  (2ttNF  t)  +  J.(NAF/f  )  cos(2tt(NF  t  +  f  )t) 
o  o  max  v  o  1  v  o  v 

-J.  (NAF  /f  )  cos (2ttNF  -  f  )t)  +  .  .  .  .}  (37) 

1  max  v  o  v 

Therefore,  for  a  frequency  multiplied  signal 


o*f 


x(f  ) 

V  V 


20  log{J.(NAF  /f  )/J  (NAF  /f  )} 
e  1  max  v  o  max  v 


(38) 


One  must  be  careful  not  to  use  inappropriate  approximations  for 
Jp(X)  if  X  >  0.1. 


To  illustrate  the  effects  of  frequency  multiplication,  we  shall 
choose  an  example  from  the  GPS  system. 

Example  1:  Let  F^  =  5,115  MHz 

y  =  2  x  10'Vg 


\ 


A 

max 

f 

v 


5g 

52  Hz 


(39) 
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Then  y*  1  (52  Hz)  =  20  log  (A  yF  /2f  )  =  20  log  (4.9  x  10'4) 
'f  max  o  v 


v  1  ( 52  Hz)  =  -66dB 

If  this  signal  is  multiplied  by  N  =  308  to  get  to  1375  MHz,  we  get 


oT  1 


(52  Hz)  =  20  logf J  (NF  A  y/f  )/J  (NF  A  y/f  )} 
1  o  max '  v  o  o  max '  v 


(40) 


since 


NF  A  y/f  =  0.31 
o  max  v 

we  must  use  the  full  values  of  J  and  J,  from  the  appendix. 

o  1 

c^fv1(52  Hz)  =  20  log  (.1483/. 997)  =  -1  6dB 


(41  ) 


(42) 


The  sideband  is  thus  only  16dB  below  the  carrier. 


Example  2: 

If  we  lower  the  vibration  frequency  to  10  Hz,  we  get 

o?f  v1(10  Hz)  =  20  logfJjd.SD/J  (1.51)}  =  +0, 87dB  (43) 

The  plus  sign  denotes  that  the  sideband  is  actually  larger  than 
the  carrier  I 

Example  3: 


If  f  =  6.55  Hz  with  A  still  5g  and  N  =  308. 
v  max 

v1( 6.55  Hz)  =  20  log{J1(2.405)/Jo(2.405) 


(44) 


Since  J  (2.405)  =  0,  1(6.55  Hz)  =  “,  i.e.,  the  carrier  disappears!! 

o  v 

Figure  5,  an  experimental  demonstration,  is  a  comparison  of  the 
power  spectrum  of  a  resonator  being  vibrated  with  its  output  frequency 
multiplied  by  a  factor  of  10,  and  the  same  resonator  with  an  unmult ipl it'd 
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output.  Since 


o?f  1 


v  (fv)(N  =  1)  =  20  log(AmaxFoY/2fv) 


(45) 


and 


X(f  ,)(N  =  10)  =  20  log(10  A  F  y/2f  ) 


max  o 


(46) 


:(N  =  10)  ~y  \n  =  1)  =  20  log(10  A  F  y/2f  )/(A  F  y/2f  ) 
v  v  °  max  o  v  max  o  v 

v1(N  =  10)  ~^fv1(N  =  1)  =  20  log(10)  =  20dB 

in  Figure  5  (N  =  10)  =  -42dB  and  (N  =  1)  =  -62dB  so 
y  L(N  =  10)  -  £  A(N  =  1)  -  20dB.  It  can  be  seen  that  as  one  increases 
tKe  multiplicative  factor,  N,  or  decreases  the  vibration  frequency,  f  , 
the  modulation  factor,  (NAF  /f  ),  will  increase.  The  oscillatory  V 
behavior  of  the  Bessel  Func?f onsVleads  to  the  growing  importance  of  higher 
order  sidebands. 

Example  4:  As  before,  if 


F  =  5.115  MHz 
o 

y  =  2  x  10'Vg 


A  =  5g 
max 

N  =  308 


(47) 


but 


f  =  52  Hz 
v 


(52) 


-lbdB 


2 (52)  =  -38dB 


The  second  order  sideband  is  only  20dB  below  the  first  order  side¬ 
band.  Before  multiplication  the  second  order  sideband  was  138dB  below 
the  carrier  and  72dB  below  the  first  order  sideband.  When  fv  =  6.55  Hz, 
the  carrier  disappears  and  the  second  order  sideband  is  only  1.6dB  below 
the  first  order  sideband. 

It  should  be  noted  that  it  is  not  the  introduction  of  the  multi¬ 
plicative  factor  which  causes  the  redistribution  of  the  power  from  the 
carrier  to  the  sidebands,  but  rather  the  increase  in  absolute  frequency. 

is  a  function  of  (NAF  /f  )  which  can  be  written  as  (yA  NF  /f  ). 

max  «  ”■ —  ~  •• 


max 
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If  NF  is  relabeled  as  F,  is  then  a  function  of  (yA  F/f  ).  It 

does  not  matter  whether  F  is  obtained  by  F  =  FQ  directly3and  ii  =  1  or 
by  F0  being  multiplied  up  to  F.  As  a  matter  of  fact,  a  fractional  value 
for  N  exhibits  the  same  vibration  response.  The  only  wav  to  reduce  the 
absolute  frequency  excursion,  I.e.,  improve  the  vibration  insensitivity, 
is  to  reduce  the  factor  y. 


Example  5: 

A  vibration  resistant  oscillator  is  needed  at  1575  MHz,  The 
designer  has  the  option  of  using  a  5.115  MHz  bulk  wave  crystal  oscillator 
multiplied  up  to  1575  MHz  or  using  a  1575  MHz  surface  acoustic  wave  (SAW) 
oscillator.  Both  oscillators  exhibit  an  acceleration  sensitivity  of 
2  x  10  9/g.  Which  oscillator  should  the  designer  choose? 

In  example  1  we  calculated  that  NF0Amaxy/fv  =  0.31  which  implied 
that  V1 (52  Hz)  =  -16dB. 

For  the  SAW  oscillator,  at  the  same  vibration  frequency  FoNiax^v 
is  agair  0.31,  which  implies  again  that  ^(52  Hz)  is  again  -16dB. 

Thus  we  see  that  the  vibration  induced  sidebands  are  identical  for 
the  two  cases.  Therefore,  the  designer  should  base  his  decision  on  other 
factors,  such  as  long  term  stability,  or  else  choose  the  device  that  is 
available  with  a  lower  vibration  sensitivity. 

The  result  of  example  6  is  easier  to  understand  in  the  time  domain, 
e.g.,  as  shown  in  equation  5  and  Figure  1.  The  frequency  of  the  5.115  MHz 
oscillator  will  vary  by  YFoA^x  =  0.05  Hz. 

Upon  multiplication  by  308,  both  FQ  and  Fmax  increase  by  a  factor 
of  308,  to  1575  MHz  and  15.75  Hz,  respectively.  The  SAW  oscillator 
frequency  will  vary  about  1575  MHz  by  yFo^max  =  -*-5.75  Hz,  which  is  the 
same  as  the  variation  in  the  multiplied  bulk  wave  oscillator  frequency. 

COMPLEX  VIBRATION 


If  the  resonator  is  being  vibrated  by  a  superposition  of  simple 
harmonic  vibrations  the  situation  may  become  more  complicated.  The 
frequency  of  the  resonator  becomes,  in  general: 


o  nmax 

n=  1 


F(t)  =  F  +  y  AF  cos  (  2tt  f  t) 
and  the  phase  is  therefore, 


$(t)  =  2ttF  t  +  l  (AF  /  f  )  sin(2rf  t), 
o  max  vn'  vn 

n=i 


(48) 


(49) 
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therefore , 


V(t)  =  V  {cos  (2ttF  t)  +  I  (AF  /f  )  sin(2iTf  t)} 
o  o  n=2  max  v  vn 


(30; 


which  becomes 


V(t)  =  V  (cos(2ttF  t)  cos(  l  (AF  /f  )  sin(2irf  t)) 
o  o  n=1  max  v  vn 


-sin(27;F  t)  sin(  I  (AF  /f  )  sin(2irf  t)} 


n=l  • 


max  v 


vn 


(51) 


If,  for  simplicity,  we  choose  only  2f  's,  we  have 


V(t)  =  V  {cos(2ttF  t)  cos(AF,  /f  ,)  sin(2TTf  t)  +  (AF„ 

o  o  lmax  vl  v  2max 

-  sin(2TrFot)  sin(AFlmax/fvl)  8in(2irfvlt)  +  (AF2/fv2) 


/fv2)  sin(2TTfy2t)) 
sin(2iTfv21:)) )  (52) 


This  becomes  unwieldly  very  quickly  unless  we  make  some  simplifying 
assumptions.  If  we  assume  for  the  moment  that  AFp,ax/f  is  small  com¬ 
pared  to  1  we  may  drop  all  terms  with  Bessel  functions  of  order  ^2  and 
all  terms  in  J^2.  V(t)  becomes  (see  Appendix  III) 


V(t)  =  V{J  (AF.  /f  .)  J  (AF  / f  ,) 
o  lmax  vl  o  2max  v2 


■J.  (AF,  /f  . )  J  (AF„  /f  „)  (cos  2  (ttF  -f 
1  lmax  vl  o  2max  v2  o 


J  (AF-  /f  „)  J.  (AF  /f  .)  (cos(2tt(F  -f  „)t 
o  2max  v2  1  max  v2  o  v  2 


cos(2ttF  t) 
o 

vl)t)  -  cos(2TT(Fo+fv2)t) 

)  -  cos(2m(F  +f  )t ) 1  ■  (511 
o  v2 


For  the  approximation  we  have  made  this  reduces  to  the  linear  sum  of  all  the 
individual  sidebands.  For  situations  where  AFnmax/fV2>l  the  result  is 
much  more  complex.  In  that  case,  less  frequency  modulation  is  required 
to  mask  the  carrier. 

RANDOM  VIBRATIONS 

In  this  case  the  vibratory  power  is  distributed  over  a  range  of 
frequencies.  The  worst  case  condition  is  described  by  a  power  spectral 
density  plot  of  the  vibration  envelope. 

The  power  spectral  density  of  a  vibration  is  the  mean  square  value 
of  the  acceleration  per  unit  filter  bandwidth  (1  Hz).  The  power  spectral 
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density  is  expressed  in  units  of  g2/Hz.  The  rms  value  of  the 
acceleration  within  a  frequency  band  between  and  f2  is,  with  G  being 
the  power  spectral  density  of  acceleration: 


A* 


/  G(f)df 


(54) 


The  value  of  the  maximum  acceleration  at  a  given  vibration  frequency 
(actually  a  1  Hz  band  of  frequencies)  is 


A  (f)  =  V'T1' TT2  <7)7  (1  Hz) 

max  rms  ’ 


(55) 


The  ratio  of  the  power  in  one  of  the  1  Hz  wide  sidebands  to  the 
carrier  power,  from  the  last  section,  is  therefore  (expressed  in  dB) 

c2fvl(f)  -  10  log  Alf)i-!,Fo/I))/(Jii('rr'  *(f)rn,jro/f))J 

oTv1(f)  =  20  log^W  Arms(f)YFo/f);  1  Hz 
Example  1 :  Assume 


(56) 


F  =  5.115  MHz 
o 

Y  =  2  x  10  9/g  . 


(57) 


G ( f )  is  as  in  Figure  6. 

oT  7f)  =  20  log (7 .23  x  10_A  (f)/f)  , 

v  rms 


(58) 


and  is  as  shown  in  Figure  7. 


CONCLUSION 


The  preceding  report  describes  the  effects  of  vibration  on  crystal 
resonators.  The  2  x  10  9per  g  used  in  the  examples  is  the  state-of-the-art 
for  commercially  available  oscillators  in  1979.  Techniaues  in  development 
such  as  the  use  of  doubly  rotated  cuts8’9,  compound  vibrators10,  and 
improved  mounting  and  bonding  techniques11  promise  to  provide  orders  of 
magnitude  improvement  in  vibration-resistance  in  the  future. 
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APPENDIX  I 


FOURIER  DECOMPOSITION  OF  cos  (8  sin  x)  and  sin(£sin  x) 

Problem:  Fourier  decompose  Cos^(8sin  x)  and  sin(8sin  x) .  Consider 

the  Maclaurin  series  expansion  of  e  and  e-^/2t> 

OO 

eSt/2  =  l  (1/r! ) (St/2) r  (Al) 

r=o 


and 


e~6/2c  =  I  (l/sI)(-B/2t)s  t  +  0 
s=o 


(A2) 


If  we  multiply  (Al)  by  (A2)  we  have 
3t/2.  -B/2t)  _  '  (3/2) (t-t-1 )  _  - 


DO  DO 


(ew  ~)(e 


=  e 


l  ((-l)S/(r:s:))(3/2)r+str"s(A3) 


r-o  s=o 


ihe  coefficient  of  the  t  term  is 


1  -  32 / 2 2  +  3“/(22-42)  -6s/(22-42*62)  +  .  .  .  -  J  (S) 


(A4) 


nth 

The  coefficient  of  the  t  term  is 


(3/2)  /n!  {1  -  (3/2) 2/  (1* (n+1) )  +  (3/2)“/ (1*2- (n+1) (n+2))>  . 


=  Jfl(3)  (A5) 


tneretore 


e2/2(t-t  l)  =  j  (8)  +  tj  (3)  +  t2 J9 (3)  +  t_1J  .(3) 
old  -1 


+  t~2j  (3)  +  .  .  .  =  y  tnj  (3) 

-  n— »  n 


(A6) 


lx 


If  we  substitute  e  ‘  for  t  we  get 


ix  -ix.  _  (i  3sin  x) 


(3/2) (t-t~! )_  (3/2)  (e  -e~1X)  =  el 


=  cos(3sin  x)  +  i  sin(3sin  x) 


(A7) 
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so 


cos(3sin  x)  +  i  sin(gsin  x)  =  J  einXJ  (g) 


=  l  {elnXJ  (g)  -  e  lnXJ  (g)}  +  J  (;-)  <  Ah) 

n=-  *  n 


Using  the  identity 


J_n(B)  =  (-1)  Jn(B) 


we  get 


cos(gsin  x)  +  i  sin(gsin  x)  =  J  (g)  +  £  (elnX+(-l)ne  1§  X).I  (g) 

°  n-1  n 


Using  the  relationship 


sin  x  =  (elx  -  e  lx)/2i 


,  ix  .  -ix. 
cos  x  =  (e  +  e  )/2 


we  get  from  equation  (A10) 


cos(gsin  x)  !  i  sin(gsin  x)  =  J  (6)  +2  J  J  (3)  cos(2kx) 

°  k=l 

OO 

+  2i  l  ^2k+l sin( (2k+l)x) 
k=o 

If  we  separate  real  and  imaginary  parts  we  have  our  result 


cos(gsin  x)  =  J  (g)  +2  [  J„,  (3)  cos(2kx) 

0  k=l 


sin(8sin  x)  =  2  [  J2k+1^  sin((2k+l)x) 


(A  1 r?) 
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APPENDIX  II 


BESSEL  FUNCTIONS 
Bessel  Functions  are  defined  as 

CO 

j  (B)  =  I  ((-i)k(B/2)2k+p/kir(p+k+i)) 

P  k=o 

if  p  is  an  integer,  the  gamma  function  is  then 

r(p+kl)  =  (p+k)! 

therefore,  for  irtteger  p 

OO 

J  (B)  =  l  ( C-l)kC6/2)2k+p/k! (P+k) : ) 

P  k=o 

Useful  relations  are 

J  (0)  =  1  and  J  (B<0.1)  =  1 
o  o 


(Bl) 


(B2) 


(B3) 


(B4) 


and 


J  x1(0)  =  0  and  J. (8<0.1)  ~  6/2  (B5) 

P>1  1 

Figure  8  is  a  plot  of  J0(B),  J, (B)  and  J-CB)  for  B  =  0  to  10. 

Figure  9  is  a  plot  of  J ( 6 ) / J_  ( 6 )  and  Jo  (S) / J0  (3)  for  B  =  0  to  0.1. 

It  can  be  seen  that  Ji(8)/J0(B)  is  -  B/2  and  that  J 2 ( B ) / JQ ( 6 )  is  much 
less  than  (B/2)  in  that  region. 

J  (B)  has  its  first  zero  at  2.405.  This  is  where  J-i  (8)/Jo(8)  goes 

.  oQ  .  ■L 

to  infinity. 

Table  I  is  a  listing  of  values  for  J  ,  J..  ,  J,,  J./J  and  J~/J  for 
3  =  0.001  to  2.  °  °  ° 


22 


J 


FIGURE  8  -  BESSEL  FUNCTIONS  OF  ORDER  0.  1.  AND 


— 0.  05  i.  PI GURE  9  -  J i  OO  / <X)  and  J 2  ^  ^'10 


TABLE  1  -  VALUES  OF  BESSEL  FUNCTIONS 


X 

JO  (X) 

Jl(X) 

J2(X) 

Jl/JO 

J2/J0 

O.Ol 

1.0000 

0.0050 

0.0000 

5.000E-03 

1.250E-05 

0.02 

0.9999 

0.0100 

0.0000 

1.000E-02 

5.000E-05 

0.03 

0.9993 

0.0150 

0.0001 

1. 500E-02 

1. 125E-04 

0.04 

0.9996 

0.0200 

0.0002 

2  .000E-02 

2  .001F-04 

0.05 

0.9994 

0.0250 

0.0003 

2 .501E-02 

3. 126E-04 

0.06 

0.9991 

0.0300 

0.0004 

3 . 001E-02 

4 . 503E-04 

0.07 

0.9988 

0.0350 

0.0006 

3.502E-02 

6.130E-04 

0.08 

0.9984 

0.0400 

0.0008 

4 . 003E-02 

3 .009E-04 

0.09 

0.9980 

0.0450 

0.0010 

4.505E-02 

1.014E-03 

0.10 

0.9975 

0.0499 

0.0012 

5. 006 E- 02 

1.252E-03 

0.11 

0.9970 

0.0549 

0.0015 

5 . SG8E-02 

1. 516E-03 

0.12 

0.9964 

0.0599 

0.0018 

6 .011E-02 

1.804E-03 

0.13 

0. 9958 

0.0649 

u. 0021 

6 .514E-02 

2  .118E-C2 

0.14 

0.9951 

0.0698 

0.0024 

7.C17F-02 

2.458E-03 

0.15 

0.994-4 

0.0748 

0.0028 

7 . 521E-02 

2.323E-03 

0.16 

0.9936 

0.0797 

0.0032 

8.026E-02 

3. 214E-0  3 

0.17 

0.9928 

0.0847 

0.0036 

8. 531c- 02 

3.630E-03 

0.16 

0.9919 

0.0896 

0.0040 

9  .0 37E-9 2 

4 .072E-03 

0.19 

0. 9910 

0.0946 

0.0045 

9.543E-J2 

4 . 5 4  OE-O  3 

0.20 

0.9900 

0.0995 

0.0050 

1 . 005E-C 1 

5 .034E-03 

0.21 

0.9890 

0.1044 

0.0055 

1.056E-01 

5 . 553E-0  3 

0.22 

0.9879 

0.1093 

0.0060 

1.107E-01 

6 . C99E-C  3 

0.23 

0.9868 

0.1142 

0.0066 

1. 158E-01 

6 .671E-03 

0.24 

0.9857 

0.1191 

0.0072 

1. 209E-01 

7 . 27  0E-C  j 

0.25 

0.9844 

0.1240 

0.0078 

1.260E-01 

7.BS5E-03 

0.2  6 

0.9832 

0.1289 

0.0084 

1.311E-01 

8 . 5  4  6  E-  0  3 

0.27 

0.9819 

0.1338 

0.0091 

1 . 362E-01 

9.225E-C3 

0.28 

0.9805 

0.1386 

0.0097 

1.414E-01 

S.  930E-0  3 

0.29 

0.9791 

0.1435 

0.0104 

1.465E-01 

1 .066E-02 

0.30 

0.9776 

0.1483 

0.0112 

1.517E-01 

1 . 142E-02 

0.31 

0.9761 

0.1531  • 

0.0119 

1. 569E-01 

1. 221E-02 

0.32 

0.9746 

0.1580 

0.0127 

1.621E-01 

1 . 302E-02 

0.33 

0.9730 

0.1628 

0.0135 

1.673E-01 

1 . 366E-C2 

0.34 

0.9713 

0.1676 

0.014  3 

1.725E-01 

1 .473E-02 

0.35 

0.9696 

0.1723 

0.0152 

1.777E-01 

1  .  5  C3F-02 

0.36 

0.9679 

0.1771 

0.0160 

1.830E-01 

1 .6  56E-02 

0.37 

0.9661 

0.1819 

0.016  9 

1. 882E-01 

1  .  751C-02 

0.38 

0.9642 

0.1806 

0.0178 

1.935E-01 

1. 850E-02 

0.39 

0.9623 

0.1913 

0.0188 

1.988E-01 

1.951 E-C  2 

0.40 

0.9604 

0.1960 

0.0197 

2 . 04 1E-0 1 

2  .0 55E-02 

0.41 

0.9584 

0.2007 

0.0207 

2  .094  F.-0 1 

2 . 162E-0  2 

0.42 

0.9564 

0.2054 

0.0217 

2 .148E-01 

2 . 272E-02 

0.43 

0.9543 

0,2101 

0.0228 

2 .201E-01 

2 . 38  5 E- 02 

0.44 

0.9522 

0.2147 

0.023B 

2.255E-01 

2 .501 E-02 

0.45 

0.9500 

0.2194 

0.0249 

2 .309E-01 

2 .620E-02 

0.46 

0.9478 

0.2240 

0.0260 

2 .3G3E-01 

2 .74  2F-02 

0.47 

0.9455 

0.2286 

0.0271 

2.417F  -01 

2 .8  67E-C2 

0.48 

0.9432 

0.2332 

0.0283 

2.472E-01 

2 . 995E-02 

0.49 

0.9409 

0.2377 

0.0294 

2 .527E-01 

3. 127E-02 

0.50 

0.9385 

0.2423 

0.0306 

2 . 582E-01 

3. 261E-02 
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TABLE  1  (CONTINUED) 


X 

JO(X) 

J1(X) 

J2(X) 

Jl/JO 

J2/J0 

0.51 

0. 9360 

0.2468 

0.0318 

2  .637E-01 

3. 399E-02 

0.52 

0.9335 

0.2513 

0.0330 

2.692E-01 

3.540E-02 

0.53 

0.9310 

0.2558 

0.0343 

2.743E-01 

3.684E-02 

0.54 

0.9284 

0. 2603 

0.0356 

2 .803E-01 

3. 831 E- 02 

0.  55 

0.9258 

0.2647 

0.0369 

2. 860 E- 01 

3. 982E-02 

0.56 

0.9231 

0.2692 

0.0382 

2.916E-01 

4.137E-02 

0.57 

0.9204 

0.2736 

0.0395 

2.972E-01 

4.294E-02 

0.58 

0.9177 

0.2780 

0.0409 

3 . 029E-01 

4.4  55E-02 

0.59 

0.9149 

0.2823 

0.0423 

3.086E-01 

4 . 620  £-02 

0.60 

0.9120 

0.2867 

0.0437 

3.144E-01 

4.788E-02 

0.61 

0.9091 

0.2910 

0.0451 

3. 201E-01 

4.959E-02 

0.62 

0.9062 

0.2953 

0.0465 

3. 259E-01 

5. 135 E- 02 

0.6  3 

0.9032 

0.2996 

C.0480 

3. 317E-01 

5.31 3E-02 

0.64 

0.9002 

0.3039 

0.0495 

3.376E-01 

5.4  96E-02 

0.65 

0.8971 

0.3081 

0.0510 

3.435E-01 

5 .682E-02 

0.66 

0.8940 

0.3124 

0.0525 

3.494E-01 

5 . 872E-02 

0.67 

0.8909 

0.3166 

0.0540 

3. 553E-01 

6.066E-02 

0.68 

0.8877 

0.3207 

0.0556 

3.613E-01 

6.264E-02 

0.69 

0.8845 

0. 3249 

0.0572 

3.673E-01 

6.466E-02 

0.  70 

0.8812 

0.3290 

0.0588 

3.733E-01 

6 .671E-02 

0.71 

0.8779 

0.3331 

0.0604 

3.794E-01 

6 . 88  IE-02 

0.72 

0.8745 

0.3372 

0.0620 

3.855E-01 

7. 095 E- 02 

0.73 

0.8711 

0.3412 

0.0637 

3. 917 E- 01 

7.313E-02 

0.74 

0.8677 

0.3452 

0.0654 

3.979E-01 

7.535E-02 

0.75 

0.8642 

0.3492 

0.0671 

4.041E-01 

7. 761 E- 02 

0.76 

0.8607 

0.3532 

0.0688 

4.104E-01 

7.992E-02 

0.77 

0.8572 

0.3572 

0.0705 

4.167E-01 

8 . 227E-02 

0.78 

0.8536 

0.3611 

0.0723 

4.230E-01 

8.466E-02 

0.79 

0.8500 

0.3650 

0.0740 

4.294E-01 

8. 710 E- 02 

0.80 

0.8463 

0.3688 

0.0758 

4.3  58E-01 

8. 959E-02 

0.81 

0.8426 

0.3727 

0.0776 

4 . 4  23E-01 

9.212E-02 

0.82 

0.8388 

0.3765 

0.0794 

4.488E-01 

9. 470 E- 02 

0.83 

0.8350 

0.3803 

0.0813 

4.554E-01 

9 . 733E-02 

0.84 

0.8312 

0.3840 

0.0831 

4.620E-01 

1.000E-01 

0.85 

0.8274 

0.3878 

0.0850 

4.687E-01 

1.027E-01 

0.86 

0.8235 

0. 3915 

0.0869 

4.754E-01 

1.055E-01 

0.87 

0.8195 

0.3951 

0.0888 

4.821E-01 

1.083 £-01 

0.88 

0.8156 

•  0.3988 

0.0907 

4.889E-01 

1.112E-01 

0.89 

0.8116 

0.4024 

0.0926 

4.958E-01 

1.141E-01 

0.90 

0.8075 

0.4059 

0.0946 

5.027E-01 

1.171E-01 

0.91 

0.8034 

0.4095 

0.0966 

5.097E-01 

1 . 202E-01 

0.92 

0.7993 

0.4130 

0.0985 

5. 167 E- 01 

1. 233E-01 

0.93 

0.7952 

0.4165 

0.1005 

5  .2  38E-01 

1 . 264E-01 

0.94 

0.7910 

0.4200 

0.1025 

5.309E-01 

1 . 296E-0  1 

0.95 

0.7868 

0.4234 

0.1046 

5.381E-01 

1.329E-01 

0.96 

0.7825 

0.42C8 

0. 1066 

5.454E-01 

1.362E-U1 

0.97 

0.7783 

0.4302 

0.1087 

5. 527E-01 

1.396E-01 

0.98 

0.7739 

0.4335 

0.1107 

5 .6  01E-01 

1.431E-01 

0.99 

0.7696 

0.4358 

0.1128 

5.676E-01 

1 . 466E-01 

1.00 

C.7652 

0.4401 

0.1149 

5. 751E-01 

1 . 502E-01 
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X 

JO  (X) 

J1  (X) 

J 2  (X) 

Jl/JO 

J2/J0 

1.01 

0.7608 

0.4433 

0.1170 

5.627E-01 

1.538E-01 

1.02 

0.7563 

0.4465 

0.1191 

5 .903E-01 

1.575E-01 

1.03 

0.7519 

0.4497 

0.1213 

5. 981 E- 01 

1.613E-01 

1.04 

0.7473 

0.4528 

0.1234 

6 .059E-01 

1.651E-01 

1.05 

0.7428 

0.4559 

0.1256 

6 .138E-01 

1.691E-01 

1.06 

0.7382 

0.4590 

0.1278 

6 .217E-01 

1.731E-01 

1.07 

0.7336 

0.4620 

0.1299 

6 .2S8E-01 

1.771E-01 

1.08 

0.7290 

0.4650 

0.1321 

6  .37  9E-0 1 

1 .813E-01 

1.09 

0.7243 

0.4680 

0.1343 

6 .461E-01 

1.855E-01 

1.10 

0.7196 

0.4709 

0.1366 

6 . 544E-01 

1 .898E-01 

1.11 

0.7149 

0.4738 

0.1388 

6  .6  28E-01 

1 . 941E-01 

1.12 

0.7101 

0.4767 

0.1410 

6 .712E-01 

1. 986E-01 

1.13 

0.7054 

0.4795 

0.1433 

6 . 798E-0 1 

2.031E-01 

1.14 

0.7006 

0.4823 

0.1456 

6 .884E-01 

2  .0  78E-0 1 

1.15 

0.6957 

0.4850 

0.1478 

6 .972E-01 

2.125E-01 

1.16 

0.6909 

0.4878 

0.1501 

7.060E-G1 

2 .173E-01 

1.17 

0.6860 

0.4904 

0.1524 

7.150E-01 

2.222E-01 

1.18 

0.6810 

0.4931 

0.1547 

7.240E-01 

2 .272E-01 

1.19 

0.6761 

0.4957 

0.1570 

7. 312 E- 01 

2 .32 3 E-0  1 

1.20 

0.6711 

0.4983 

0.1593 

7.425E-01 

2 . 37  4  E-0  1 

1.21 

0.6661 

0.5008 

0.1617 

7. 518E-01 

2  .4  27  E-0  1 

1.22 

0.6611 

0.5033 

0.1640 

7.613E-01 

2 .48 1E-0  1 

1.23 

0.6561 

0.5058 

0.1664 

7 . 710E-0 1 

2 .536E-01 

1.24 

0.6510 

0.5082 

0.1687 

7.807E-01 

2 . 59  2F.-0  ] 

1.25 

0.6459 

0.5106 

0.1711 

7.906  E-0 1 

2 .649E-01 

1.26 

0.6408 

0.5130 

0.1735 

8.005E-01 

2 . 707E-0  1 

1.27 

0.6356 

0.5153 

0.1758 

8. 107E-01 

2 .766E-01 

1.28 

0.6305 

0.5176 

0.1782 

8. 209E-01 

2  . 8  27  E-  0  ] 

1.29 

0.6253 

0.5198 

0.1806 

8.313D-01 

2 . 8  8  9E-0  1 

1.30 

0.620] 

0.5220 

0.1830 

8.419E-01 

2 .952E-01 

1.31 

0.6149 

0.5242 

0.1854 

8.525E-01 

3. 01CE-01 

1.32 

0.6096 

0.52C3 

0.1878 

8 . 6  34  E-0 1 

3. 081E-01 

1.33 

0.6  04  3 

0.5284 

0.1903 

8. 744E-01 

3. 148E-0  1 

1.34 

0.5990 

0.5305 

0.1927 

8.855E-C1 

3. 217E-01 

1.35 

0.5937 

0.5325 

0.1951 

8 . 968E-0  1 

3.  286  E-0  ] 

1.36 

0.5884 

0.5344 

0.1976 

9  .083E-01 

3  .  35 8 E-0  1 

1.37 

0.5830 

0.5304 

0.2000 

9  . 2  0  0  E-0 1 

3 . 4  30  E-0  1 

1.38 

0.5777 

0.5383 

0.2024 

9  .3180-01 

3. 50 5 E-0 1 

1.39 

0.5723 

0.540] 

0.2049 

9  .438E-01 

3. 581E-01 

1.40 

0.5669 

0.5419 

0.2074 

9  .561E-01 

3.658F-01 

1.41 

0.5614 

0.5437 

0.2098 

9  .6  05E-01 

3.  7  37  E-0 1 

1.42 

0.5560 

0.5455 

0.2123 

9 . 8 1 lE-0 ) 

3. 8) EE-01 

1.43 

0. 5505 

0.5472 

0.2147 

9  .9  39E-01 

3. 901E-0J 

1.44 

0.5450 

0.5488 

0.2172 

1.007E  00 

3 .9850-0  1 

1.45 

0.5395 

0.5504 

0.2197 

1.020F  00 

4 .072E-0  1 

1.46 

0.5340 

0.5520 

0.2222 

1.0341'  00 

4 .  160 E-0  1 

1.47 

0.5285 

0.5536 

0.2246 

1.047E  00 

4 .251E-01 

1.48 

0.5230 

0.5551 

0.2271 

1.061E  00 

4  .  3  4  3  E  -  0  3 

1.49 

0.5174 

0.5565 

0.2296 

1.076E  00 

4  .  4  38E-0 3 

1.50 

0.5118 

0.5579 

0.2321 

1.090C  00 

4  .  5  34  E -0  i 
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JO(X)  J1 (X)  J2  (X)  31/JO 


J2/J0 


1.51 

0.5062 

1.52 

0. 5006 

1.53 

0.4950 

1.54 

0.4894 

1.55 

0.4838 

1.56 

0.4781 

1.57 

0.4725 

1.  58 

0.4668 

1.59 

0.4611 

1.60 

0.4554 

1.61 

0.4497 

1.62 

0.4440 

1.63 

0.4383 

1.64 

0.4325 

1.65 

0.4268 

1.66 

0.4210 

1.67 

0.4153 

1.68 

0.4095 

1.69 

0.4038 

1.70 

0.3980 

1.71 

0.3922 

1.72 

0.3364 

1.73 

0.3806 

1.74 

0.3748 

1.75 

0.3690 

1.76 

0.3632 

1.77 

0.3574 

1.78 

0.3516 

1.  79 

0.3458 

l.CO 

0.3400 

1.81 

0.3342 

1.82 

0.3284 

1.83 

0.3225 

1.84 

0.3167 

1.85 

0.3109 

1.86 

0.3051 

1.87 

0.2993 

1.88 

0.2934 

1.89 

0.2870 

1.90 

0.2818 

1.91 

0.2760 

1.92 

0.2702 

1.93 

0.2644 

1.94 

0.2586 

1.95 

0.2520 

1.90 

0.2470 

1.97 

0.2412 

’  .98 

0.2354 

1.99 

0.2297 

2.00 

0.2239 

0.5593 

0.2346 

0.5607 

0.2371 

0.5620 

0.2395 

0.5632 

0.2420 

0.5644 

0.2445 

0.5656 

0.2470 

0.5667 

0.2495 

0.5678 

0.2520 

0.5689 

0.2545 

0.5699 

0.2570 

0.5709 

0.2595 

0.5718 

0.2619 

0.5727 

0.2644 

0.5735 

0.2669 

0.5743 

0.2694 

0.5751 

0.2719 

0.5758 

0.2743 

0.5765 

0.2768 

0.5772 

0.2793 

0.5778 

0.2817 

0.5783 

0.2842 

0.5788 

0.2867 

0.5793 

0.2891 

0.5798 

0.2916 

0.5802 

0.2940 

0.5805 

0.2964 

0.5808 

0.2989 

0.5811 

0.3013 

0.5813 

0.3037 

0.5815 

0.3061 

0.5G17 

0.3086 

0.5818 

0.3110 

0.5818 

0.3134 

0.5819 

0.3157 

0.5818 

0. 3181 

0.5818 

0.3205 

0.5817 

0.3229 

0.5816 

0.3252 

0.5814 

0.3276 

0.5812 

0.3299 

0.5809 

0.3323 

0.5806 

0.3346 

0.5803 

0. 3369 

0.5799 

0.3392 

0.5794 

0.3415 

0 . 57  SO 

0.3438 

0.5785 

0.3461 

0.5779 

0.3483 

0.5773 

0.3506 

0.5767 

0.3528 

1.105E  00 
1.120E  00 
1.135E  00 
1.151E  00 
1.167E  00 
1.183E  00 
1.  200  E  00 
1.216E  00 
1.234E  00 
1.251E  00 
1.269E  00 
1.288E  00 
1.307E  00 
1.326E  00 
1.346E  00 
1.366E  00 
1.387E  00 
1.408E  00 
1.429E  00 
1.452E  00 
1.475E  00 
1.498E  00 
1.522E  00 
1.547E  00 
1.572E  00 
1.598E  00 
1.625E  00 
1.653E  00 
1.681E  00 
1.710E  00 
1 . 74 IE  00 
1.772E  00 
1.804E  00 
1.837E  00 
1.872E  00 
1.  907E  00 
1.944E  00 
1.982E  00 
2.021E  00 
2.062E  00 
2.105E  00 
2.149E  00 
2.195E  00 
2.242E  00 
2.292E  00 
2.344E  00 
2.398C  00 
2.455E  00 
2.514E  00 
2.576E  00 


4 .634E-01 
4. 735E-01 
4.839E-01 
4 . 946E-01 
5.055E-01 
5.166E-01 

5.  281E-01 
5. 399E-01 
5.519E-01 
5.643E-01 
5.770E-01 
5.900E-01 
6.033E-01 
6.171E-01 
6.312E-01 
6.457E-01 
6.606E-01 

6 .  759E-0 1 
6.917E-01 
7.079E-01 
7. 246 E- 01 
7.418E-01 
7.596E-01 
7  .778E-01 
7 .967 E- 01 
8.161E-D1 
8.362E-01 
8.569E-01 
8.783E-01 
9.005E-01 
9.233E-01 
9.470E-01 
9.715E-01 
9.969E-01 
1.023E  00 
1.051E  00 
1.079E  00 
1.108E  00 
1.139E  00 
1.171E  00 
1.204E  00 
1.238C  00 
1.274E  00 
1.312E  00 
1.351E  00 
1.392E  00 
1.435E  00 
1.480E  CO 
1.527E  00 
1.576E  00 
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I 


X 

JO  (X) 

J1  (X) 

2  .01 

0  .2181 

0.5761 

2  .02 

0.2124 

0.5754 

2.03 

0.2066 

0.5746 

2.04 

0.2009 

0.5738 

2  .05 

0.1951 

0.5730 

2.06 

0.1894 

0.5721 

2.07 

0.1837 

0.5712 

2  .08 

0.1780 

0.5703 

2.09 

0.1723 

0.5693 

2.10 

0.1666 

0.5683 

2.11 

0.1609 

0.5672 

2.12 

0.1553 

0.5661 

2.1: 

0.1496 

0.5650 

2.14 

0.1440 

0.5638 

2  .15 

0.1383 

0.5626 

2.16 

0.1327 

0.5614 

2.17 

0.1271 

0.5601 

2.18 

0.1215 

0.5587 

2.19 

0.1159 

c.5574 

2.20 

0.1104 

0.5560 

2.21 

0.1048 

0.5545 

2.22 

0.0993 

0.5530 

2.23 

0.0937 

0.5515 

2.24 

0.0882 

0.5500 

2.25 

0.0827 

0.5484 

2.26 

0.0773 

0.5468 

2.27 

0.0718 

0.5451 

2  .28 

0.0664 

0.5434 

2.29 

0.0609 

0.5416 

2.30 

0.0555 

0.5399 

2.31 

0.0502 

0.5381 

2.32 

0.0448 

0.5362 

2.33 

0.0394 

0.5343 

2.34 

0.0341 

0.5324 

2.35 

0.0288 

0.5305 

7.36 

0.0235 

0.5285 

2.37 

0.0182 

0.5265 

2.38 

0.0130 

0.5244 

2.39 

0.0077 

0.5223 

2.40 

0.0025 

0.5202 

2.41 

-0.0027 

0.5180 

2.42 

-0.0079 

0.5158 

2.43 

-0.0130 

0.5136 

2.44 

-0.0181 

0.5113 

2.45 

-0.0232 

0.5091 

2.46 

-0.0283 

0.5067 

2.47 

-0.0334 

0.5044 

2.48 

-0.0384 

0 . 5020 

2.49 

-0.0434 

0.4996 

2.  50 

-0.0484 

0.4971 

J2  (X) 

Jl/JO 

J  2/JO 

0.3551 

2.6  411  00 

1.6286 

00 

0.3573 

2 .769!'  00 

1 . 6  8  2  6 

co 

0.3595 

2.7811  C'j 

1  .  74  0  6 

0  0 

0.3617 

2  .857L  CO 

1.8016 

CO 

0.3639 

2  .9  36  f  05 

1  .  8  6  56. 

OC 

0.3661 

3.0216  00 

1.9336 

CO 

0.3682 

3.1106  00 

2  .004  6. 

cc 

0.3704 

3.2046  00 

2.081  6 

CO 

0.3725 

3.  3C4  6  00 

2.1626 

00 

0.37^6 

3.41U  CO 

2  . 2  4  9  E 

OC 

0.3767 

3.5256  CO 

2.3416 

6  0 

0.3788 

3.6466  OC 

2  . 4  4  c  r: 

cc 

0.3809 

3.7776  0  0 

2  .  5  4  C  ! 

6  l 

0.3830 

3.91  II  60 

2 .  <  c  or 

0.3850 

4.067f  CC 

2  .  7  k  3  1 

rr 

0.3871 

4.2306  CO 

2  . 5  ;  7 1 

;  -  ( 

0. 3891 

4.4066  CO 

?.o(.ll 

c  c 

0.3911 

4.5986  00 

3.  21 51 

z  c 

0.3931 

4.8086  00 

3.391! 

V. 

0.3951 

5.0386  OC 

3. 56.06 

L  i 

0.3970 

5.2916  00 

3  .  7  3*'! 

0.3990 

5.5716  00 

4. Cl?! 

f  ( 

0.4009 

5.8836  CO 

4 . 2766 

0.4028 

6.2  316  CO 

4.565! 

>>  < 

0.4  04  7 

6  .6  276  0  0 

4.6911 

0.4066 

7.0756  Uc 

5. 261  1 

Cl) 

0.4084 

7  .  5906  00 

5.6  871. 

u> 

0.4103 

8.1876  GO 

6.1611 

r  z 

0.4121 

8.8876  00 

6  .762  6 

V  ( 

0.4139 

9 .7206  CO 

7  .  4  0  36 

0.4157 

1.0736  Cl 

8 .269  6 

0.4175 

1.1976  01 

9  .  32  ;; 

’ 

0.4192 

1 .3556  0 1 

1  .  U 

0.4210 

1.5626  01 

i .  2  < r ; 

V  } 

0.4227 

1.6436  01 

1  .  4  i  ! 

0.4244 

2  .2506  0  1 

1.6.7! 

0.4261 

2.8916  M 

2  .  3  •*.  C  : 

0.4277 

4.0486  01 

_  i 

0.4294 

6.7666  01 

0  .  06.  ! 

0.4310 

-  2.0746  0> 

1  .  7  1  "  ! 

0.4326 

-  1  .  936  6  6. 

-1  .(  1.6 

i  : 

0.4342 

-6.5696  Cl 

-  5  . 

0.4357 

-3.9516  01 

-  '  .  .15?! 

,  ] 

0.4373 

-2.8216  01 

-2.41  11 

0.4388 

-2.1926  01 

-1.8!  '  ■ ! 

l  1 

0.4403 

-1  .  7  90  6  0 i 

-  1  ?  ? 

0.4418 

-  1 . 512'-  0! 

-1.  3  2  ■■ ! 

i  ■ 

0.4432 

-  1  .  30  7  1,  n 

-1.15 4 ! 

0.4446 

-1.1516  01 

- 1 .  r :  ■  .• 

0.4461 

-1.027E  01 

-9.219: 

l1  i 

TABLE  1  (CONTINUED) 


X 

JO  (X) 

J1(X) 

J2  (X) 

Jl'/JO 

2.51 

-0.0533 

0.4946 

0.4475 

-9  . 272E 

00 

2.52 

-0.0583 

0.4921 

0.4488 

-8.444E 

00 

2  .53 

-0.0632 

0.4895 

0.4502 

-7. 748E 

00 

2  .54 

-0.068  1 

0.4870 

0.4515 

-7.154E 

00 

2.55 

-0.0729 

0.4843 

0.4528 

-6.642E 

00 

2.56 

-0.0778 

0.4817 

0.4541 

-6.195E 

00 

2 . 57 

-0.0826 

0.4790 

C.4553 

-5.802E 

00 

2.58 

-0.0873 

0.4763 

0. 4  566 

-5.454E 

00 

2 . 59 

-0.0921 

0.4736 

0.4578 

-5.143E 

00 

2  .60 

-0.0968 

0.4708 

0.4590 

-4.864E 

00 

2.61 

-0.1C15 

0.4680 

0.4601 

-4.611E 

00 

2.62 

-0.1C62 

0.4652 

0.4613 

-4.3G2E 

00 

2.63 

-0.1108 

0.4624 

0. 4624 

-4.173E 

00 

2  .64 

-0.1154 

0.4595 

0.4635 

-3 . 98  IE 

00 

2.65 

-0.1200 

0.4566 

0.4646 

-3.805E 

00 

2  .66 

-0.1245 

0.4536 

0.4656 

-3.642E 

00 

2.6  7 

-0.1291 

0.4507 

0.4666 

-3. 492E 

00 

2  .68 

-0.1336 

0.4477 

0.4676 

-3.352F 

00 

2.69 

-0.1380 

0.4446 

0.4686 

-3. 222E 

00 

2 . 70 

-0.1424 

0.4416 

0.4696 

-3. 10CE 

00 

2.71 

-0 .1469 

0.4385 

0.4705 

-2.9C6E 

00 

2 . 72 

-0. 151  2 

0.4354 

0.4714 

-2.879E 

00 

2.73 

-0.1556 

0.4323 

0.4723 

-2.779E 

00 

2.74 

-0.1599 

0.4291 

0.4731 

-2  .6  8  4  E 

00 

2 . 75 

-0.1641 

0.4260 

0.4739 

-2.595E 

00 

2.76 

-0.1684 

0.4228 

0.4747 

-2  .SHE 

00 

2  .77 

-0.172b 

0.4195 

0.4755 

-2.431E 

00 

2. 78 

-0.1768 

0.4163 

0.4763 

-2 . 35 5 E 

00 

2 . 79 

-0.1809 

0.4130 

0.4770 

-2  .28  3E 

00 

2.80 

-0. 1C50 

0.4097 

0.4777 

-2 .2 14  E 

00 

2.81 

-0.1891 

0.4064' 

0.4784 

-2 . 1 4  9E 

00 

2 . 82 

-0.1932 

0.4030 

0.4790 

-2  . 0  8  6E 

00 

2.83 

-0.1972 

0.3997 

0.4796 

-2  .027E 

00 

2  .84 

-0.2012 

0.3963 

0.48C? 

-1.  970  E 

00 

2 .85 

-0.2051 

0.3928 

0.4808 

-1.915E 

00 

2  .  86 

-0.2090 

0.3894 

0.4813 

-1 .863E 

00 

2.87 

-0.2129 

0.3859 

0.4818 

-1.813E 

00 

2.88 

-0.2167 

0.3825 

0.4823 

-1  .  765E 

00 

2.89 

-0 .2  20  5 

0.3790 

0.4828 

-1. 718E 

00 

2  .90 

-0.2243 

0.3754 

0.4832 

-1.674E 

00 

2 . 91 

-0.2  28  0 

0. 371 9 

0.  4836 

-1. 63  IE 

00 

2 . 92 

-0.2317 

0. 3683 

0.4840 

-1.589E 

00 

2.93 

-0.2354 

0.3647 

0.4844 

-1 . 549E 

00 

2 . 94 

-0.2  590 

0. 36  11 

0.4847 

-1 . 51  IE 

00 

2 . 95 

-0.2426 

0.3575 

0.4850 

-1  .  473E 

00 

2  .  96 

-0.2462 

0. 3538 

0.4C5J 

-1  .  4.97E 

00 

2  .  97 

-0.2497 

0.3502 

0.4855 

-1.4C2E 

00 

2 . 98 

-0.2532 

0. 3465 

0.4857 

-1  .  3  68  r 

00 

2 . 99 

-0.2566 

0. 3428 

0.4859 

-1  .  3  36  E 

00 

3.00 

-0.2001 

0.3391 

0.4861 

-1.  304  E 

00 

J2/J0 

-8.388E  00 
-7.702E  00 
-7.125E  00 
-6.633E  00 
-6.209E  00 
-5.840E  00 
-5.515E  OC 
-5. 228E  00 
-4.971E  00 
-4.741E  00 
-4.533E  00 
-4.345E  00 
-4.173E  00 
-4.016E  00 
-3.872E  00 
-3 . 739E  00 
-3.616E  00 
-3.501E  00 
-3. 395E  00 
-3.296E  00 
-3.204E  00 
-3.117E  00 
-3.036E  00 
-2  .959E  00 
-2  -887E  00 
-2.819E  CO 
-2  .755E  00 
-2.694E  00 
-2.636E  00 
-2.532E  00 
-2  .529E  00 
-2.480E  00 
-2  .432E  00 
-2 . 387E  00 
-2.344E  00 
-2.303E  00 
-2.263E  00 
-2  .2250  00 
-2 . 1 8 9E  00 
-2.154E  00 
-2.121E  DO 
-2.089E  00 
-2  .0 58E  00 
-2.028E  00 
-1.999E  00 
-1.971E  00 
-1.944E  00 
-1.91BE  00 
-1.893E  00 
-1.869E  00 
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APPENDIX  III 


COMPLEX  VIBRATIONS 

Calculations  Leading  to  the  Result  for  Complex  Vibrations 
We  have  from  equation  52 

V(t)  =  V  {cos  (2ttF  t)  cos  ( (AF  /  f  J  sin(27rf  ,t)) 
o  o  1  vl  vi 

+  (AF2/fv2>  sin(2TTfv2t))  - 

sin(2nF  t)  sin((AF,/f  .)  sin(2TTf  t)) 
o  1  vl  vl 

+  (AF  /f  )  sin(2’if  0t)  • 

2  vi  v2 

Using  equation  (12)  on  (Cl)  and  substituting 
Xx  =  (AF  /fvl)  sin(2trft)t) 

X2  =  (AF2/fv2)  sin(2^f  2t) 

we  get 

V(t)  =  V  {cos(2ttF  t)  cosX,  cosX.-cos  (2^F  t)  sinX.  sinX 
o  o  12  o  12 

-  sin(2’F  t)  sinX,  cosX_-sin(2”F  t)  cosX,  sinX.  . 

o  12  o  12 

If  AF  / f  .  and  AF_/f  _  and  substituting  8,  for  M',/f  .  and 

,  _  .,1  vl  _  2  v2  1  1  vl 

AF 2/ f v2  we  8et 

cos (X . )  =  J  (2.)  (i  =  1,2) 

1  01 

and 

sin(X . )  =  2J , ( 2 . )  sin(2"f  ,t)(i  =  1,2)  . 

1  1  1  v  1 

Using  (C5)  and  (06)  in  (04)  we  get 


V(t)  =  V  {J  (8.)  J  (8,)  cos(2ttF  t) 
o  o  1  o  2  o 


-  4  J  (8  )  J  (3  )  sin(2irf  ,t)  sin(27Tf  ,t)  cos(2ttF  t) 
1112  vl  v2  o 


-  2  J  (6  )  J,  (3_)  sin(2TTf  _t)  sin(2TrF  t) 
1112  vl  o 


+  2  J  (6.)  2,(8,)  sin(2wf  t)  sin(2iTF  t)} 
o  1  1  2  v2  o 


The  second  term  can  be  decomposed  by 


cos  A  sin  B  sin  C  =  cos  A (-*5  cos(B+C)  +  'i(cos  (B-C) ) 


=  H  cosA  cos(B-C)  -  %  cos  A(cos(B+C)) 


=  k  { cos (A+B+C)  +  cos (A+B-C) 


-  cos(A-B-C)  -  cos (A+B+C)} 


So, using  (C8)  in  (C7),we  get 


V(t)  =  V  J  (3.)  J  (P,)  cos (2~F  t) 
o  o  1  o  2  o 


+  J  (8  )  J  (8  )(cos(2"(F  +f  ,)t  -  cos(2TrF  -f  )t) 
1  1  o  2  o  vl  o  vl 


+  J  (8.)  J,  ('?,)  (cos  (2"  (F  +  f  ,)t  -  cos  (2”F  -f  )t) 
o  1  1  2  o  v2  o  v2 


+  terms  in  J  (3  )  .1^(8.,)-. 


The  terms  in  J,J,  have  frequencies  like  t  ,  +  f  _ ,  f  ,  -  f  „ , 
......  1  1  vl  v2  vl  v2 

t  ,  -  t  ,  and  harmonics. 

v2  vl 


f 
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